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By combining analytic and numerical methods, edge states on a finite width graphene ribbon in 
a magnetic field are studied in the framework of low-energy effective theory that takes into account 
the possibility of quantum Hall ferromagnetism (QHF) gaps and dynamically generated Dirac-like 
masses. The analysis is done for graphene ribbons with both zigzag and armchair edges. The 
characteristic features of the spectrum of the edge states in both these cases are described. In 
particular, the conditions for the existence of the gapless edge states are established. Implications 
of these results for the interpretation of recent experiments are discussed. 
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I. INTRODUCTION 

Graphene is a remarkable system with many unusualproperties that was created for the first time only a few years 
agOfi (For reviews on graphene see, for example, Refs. |2l3|lj|5|.) One of such properties is an unconventional quantum 
Hall effect (QHE). Instead of the ordinary QHE, expected in the case of 2-dimensional non-relativistic electron systems, 
an anomalous quantization is observed in graphene.^^ The observation appears to be in perfect agreement with the 
theoretical predictions^^iSiiiS stating that the QHE plateaus should occur at filling factors v = ±4(|n| + 1/2) where 
n is an integer. This anomalous quantization is a direct outcome of the relativistic-like nature of the low-energy 
quasiparticles in graphene described by a Dirac theory with an internal U{A) symmetry] -'^^'^^ 

It is the C/(4) symmetry, operating in the spin and sublattice-valley spaces, that is responsible for the 4-fold 
degeneracy of the Landau energy levels and for the overall factor 4 in the filling factors of the observed plateaus, 
V = ±4(|n| + 1/2). Strictly speaking, the C/(4) symmetry is not exact but broken down to a smaller C/(2)+ x C/(2)_ 
symmetry group by the Zeeman term. The latter symmetry operates in the sublattice-valley space and does not mix 
spin-up (s = -I-) and spin-down (s = — ) states. When the magnetic field is not too strong, a relatively small Zeeman 
term does not affect the observable QHE in a qualitative way. 

When the magnetic field becomes sufficiently strong, the QHE plateaus v ~ ±1, ±4 are observed) This 
suggests that the 4-fold degeneracy of the Landau levels is lifted. The new plateaus may be explained by one of 
the following seemingly different theoretical scenarios, (i) Quantum Hall ferromagnetism (QHF) ^^i^^i^^i^°i^^ , which is 
connected with the theory of exchange-driven spin splitting of Landau levelsi^ The QHF order parameters are densities 
of the conserved charges connected with three diagonal generators of the non-abelian subgroup SU{4) C /7(4) (the 
dynamics of a Zeeman spin splitting enhancement considered in Ref. |2^ is intimately connected with the QHF). (ii) 
The magnetic catalysis (MC) scenario2ii^i2Sii^ that is based on the phenomenon of an enhancement of the density 
of states in infrared by a magnetic field, which catalyzes electron-hole pairing (leading to excitonic condensates) 
in relativistic-like systemsi ^^'^^i'^° . This scenario invokes electron-hole pairing and excitonic condensates to produce 
dynamically generated Dirac-like masses in the low-energy theory. 

Recently, by analyzing a gap equation for the propagator of Dirac quasiparticles, it has been found in Refs. [sills^ 
that the QHF and MC order parameters necessarily coexist. As will be shown in the present paper, this feature could 
have important consequences for the dynamics of edge states in the QHE in graphene. 

The study of edge states is of general interest because such states provide a deeper insight into the quantum Hall 
effecti^ Currently there exist many studies of edge states in graphene-like systemsj "'^°'^'^i'^^'^^i'^^'^'^'^^''^^i"^°'^"'^i"'^^'^^ In 
Refs. [23II38I it was found that in the presence of a magnetic field there may exist gapless modes of such states and 
they should play an important role in charge transport of graphene near the Dirac neutral point. The gapless modes 
were shown to appear when the lowest Landau level (LLL) was split by a spin gap, or in other words, by enhanced 
Zeeman splitting leading to the = plateau. Such gapless states were absent, though, in the case of a Dirac mass 

While the presence of the gapless edge states should make graphene a so-called quantum Hall metal, their absence 
should make it an insulator li^iS^ The actual temperature dependence of the longitudinal resistivity at the 1/ ^ 
plateau in Refs. [Tslfisl is consistent with the metal type. Thus, it was argued that the origin of the v = plateau is 
connected with the enhanced spin (ferromagnetic) splitting of the LLLf^^i^ 
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The conclusion of Ref. [T5II23I regarding the origin of the 1/ = plateau does not appear to be universal however. The 
recent data from Ref. [l^ reveal a clear plateau at = 0, but the temperature dependence of the diagonal component 
of the resistivity signals a crossover to an insulating state in high fields. This does not seem to support the existence 
of gapless edge states. So, one may ask whether there is indeed a Dirac-type mass gap and no spin gap in the device 
studied in Ref. [H? 

Motivated by this question and a theoretical analysis in Refs. IsiUs^ the spectrum of edge states has been recently 
studied under the assumption that the removal of sublattice and spin degeneracies in graphcne in a strong magnetic 
field is connected with the generation of both Dirac masses and spin gaps4^ The main result of that work was 
establishing a criterion for the existence of the gapless edge states in the QHE in graphene on a half-plane with a 
zigzag or armchair edge. It was concluded that the controversy between the experimental results in Refs. [T^ and[l6l 
may reflect more rich and complicated dynamics in the QHE in graphene than those considered in the QHF and MC 
scenarios. 

In the present work, we extend the analysis in Ref. [12 to the case of a finite width graphene ribbon. The characteristic 
feature in our approach is combining analytic and numerical methods in the analysis of the edge states. This allows 
to describe the main features of the quasiparticle spectrum and, in particular, to extend the criterion obtained in 
Ref. [13 to a more complicated and interesting case of a graphene sample with two boundaries, which in turn yields a 
deeper insight in physics behind this criterion. 

The paper is organized as follows. In Sec. [Til we discuss the low-energy field-theoretical model of graphene with 
MC and QHF dynamical order parameters used in the rest of the paper. An overview of the general formalism for 
studying edge states of a finite width graphene ribbon in an external magnetic field is given in Sec. lIIII The numerical 
analysis in the case of zigzag and armchair edges is presented in Sec. IIVI and |Vl respectively. The discussion of the 
main results and their experimental implications are given in Sec. IVII 



II. QHF AND MC ORDER PARAMETERS IN GRAPHENE 



For convenience of the analysis of the edge state, the chiral representation of the Dirac matrices will be used here 
(see, for example, Ref. [1), 



75 ^ z7V7V = r3®rO=(^J ) , (2) 

where f^r' are Pauli matrices and t° is the 2x2 unit matrix. The spinor field of Dirac quasiparticles = 

(V'lf^yi: V'lf ^ B J V'lf^ B : ''Pk^a^ combines the Bloch states on two valleys (i4r+ and K^) and two sublattices {A and S), 

s is the spin index. The QHF order parameters are the spin density {^''Ps'ii) and the pseudospin density {'i/'' -f^ Pg'i^) , 
with P± — {1 ± a^)/2 being projectors on states with spin directed along {+) and opposite (— ) the magnetic field. 
These order parameters are related to the chemical potentials Hs and jig, respectively. On the other hand, the MC 
order parameter is the vacuum expectation value of the Dirac mass term {^j'^Ps'^) associated with the conventional 
Dirac mass Ag (here ^ = '^^j^). 

Recently, a unifying approach, combining and augmenting both QHF and MC mechanisms, was proposed in Refs.[3ll. 
[32I . By analyzing the gap equation with a local Coulomb interaction and using a multi-parameter variational ansatz for 
the quasiparticle propagator, it was found that (i) the new MC order parameter (^'7^7^Ps^'), related to a Dirac mass 
As that breaks time reversal symmetryi^, has to be added, and (ii) the QHF and MC order parameters necessarily 
coexist. (Let us emphasize that, in the presence of an external magnetic field, the time-reversal symmetry is broken 
and a state with the vanishing As is not protected by any symmetry.) 

More precisely, it was shown in Refs. I3lll33 that for a fixed spin, the full inverse quasiparticle propagator takes the 
following general form (in the chiral representation used in this paper): 

iG-\u, u') = [{ihdt + ^ls+ A«7')7" -VF{7T--f)- As73 + A,7375] S^u ~ u'), (3) 

where tt is the canonical momentum and the parameters fis, fJ-s, Ag, and Ag are determined from the gap equation. 
Note that the full electron chemical potentials fi± include the Zeeman energy =pZ with 

Z ~ hbB = 0.67B[T] K. (4) 
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By making use of the explicit form of the spinor. wc find the following correspondence between the four types of 



order parameters and the electron densities (the spin index is omitted): 

^ -> ('I'^^') = UK^A + nK_A + n-K+B + nK_B, (5) 

fi (*^7^*) = uk^a - nx^A + nK+B - ^^i^_B, (6) 

A ^ i^-f^j^^) = nii+A - nK_A - riK+B + nK_B, (J) 

A ^ {^7^"^) ^ nK+A + nK_A~nK+B -nK_B, (8) 



where the uk+a, nK_A, nx+B and tlk^^b are the densities of quasiparticles at specified valleys and sublattices. The 
QHF order parameters associated with ^ and Jl are the total density of electrons with a given spin and the density 
imbalance between the two valleys, respectively. Note that the MC order parameter related to the conventional Dirac 
mass A describes the density imbalance between the A and B sublattices (i.e., a charge density wave as interpreted 

m Refs.miniiiiiaiii). The value of the singlet Dirac mass A [see Eq. (7)] controls a mixed density imbalance at 

the two valleys and the two sublattices. 

In terms of symmetry, these order parameters can be divided in two groups. The order parameters {'^^ Pg'^) and 
{^7^7^Ps'^) (related to fSg and Ag) with nonequal values for s = ± break the approximate C/(4) symmetry just 
like the Zeeman term, and they are singlets under the non-abelian subgroups SU{2)s C U{2)s with the generators 
Ps <Xi (— Z7'^/2, 7'^7^/2, 7^/2) . Since these singlet order parameters break no exact symmetries of the action, they are 
not the order parameters in the strict sense. Yet, because of a relative smallness of the "bare" Zeeman energy Z and a 
significant dynamical par t in and Ag, it is appropriate to talk about approximate spontaneous symmetry breaking. 
In the model in Refs. I3lll32l . these two order parameters coexist and play a crucial role in the solution corresponding 
to the V = plateau. 

The order parameters of the other type, {'^^^^P^'^) and {^^^Pg'^) , (associated with jig and A^) are triplets under 
SU{2)s C U[2)s. Each of these two order parameters describes spontaneous SU{2)s symmetry breaking down to 
U{l)s with the generator Pg ® 7^/2. In the model in Refs. [3lll3^ . these two order parameters coexist and play an 
important role in the solution corresponding to the v = ±1 plateaus (as well as to the plateaus v = ±3 and v = ±5, 
which have not been observed yet). 

III. LANDAU LEVELS AND EDGE STATES 

In accordance with the discussion in the previous section, we assume that both ferromagnetic and mass type gaps 
coexist in general. Our goal is to find the spectrum of edge states in such a theory. 

The structure of the inverse quasiparticle propagator in Eq. Q implies that in the most general case the quadratic 
part of the effective Hamiltonian for quasiparticles of spin s takes the following form in the first quantized theory: 

H,^Ho-tis- flsl^ - A,7°73^5 + As7°j^ . (9) 

The free part of the Hamiltonian reads 

Hq = Vp [aiTTx + a2TTy) (10) 

where vp — 10^ m/s is the Fermi velocity. By definition, Ui = 7'^7' and the canonical momentum is tt^ = —ihdi+eAi/ c. 
Here the vector potential is taken in the Landau gauge: Ax = —By and Ay = 0, where B is the magnitude of a 
constant magnetic field orthogonal to the xy plane of graphene. 

The parameters /Zs, jis, Ag and Ag are determined from the gap equation. In particular, in the model in Refs. [3lll3^ . 
the 51 solution (in the nomenclature of Ref. [33 ) near the Dirac neutral point, corresponding to the v = plateau, 
has the following form: 

A±=/i±=0, Ai± = T^TA A±=±M, (11) 

where for the values of magnetic fields _B < 45 T utilized in the experiments in Refs. [islfl^fTslflGl . the dynamical 
parameters A^ M are considerably larger than Z, and M > A. The value of the spin gap in this solution is AE = 
2M + 2{Z + A). Note that it is essentially larger than the spin gap AE" = 2{A + Z) in the QHF scenario-^ While 
the problem of calculating the values of /i^, /ig, Ag and Ag is not addressed in this study, it should be clear that in 
an actual device they are determined by (i) the strength of the magnetic field, (ii) the temperature, and (iii) other 
sample-specific parameters (e.g., the mobility of carriers, the size and geometry, the type of the substrate etc.). In 
practice, we analyze the spectrum of edge states in the model described by the model Hamiltonian in Eq. ([9|). 
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In general, in a finite geometry case, the magnitude of the exchange and Hartree interactions which determine the 
values of the parameters Ag , Ag , /x^ , and jig is likely to vary with the distance from the edge and should be calculated 
in a self-consistent way. The present study of the edge states is done assuming uniform gaps and uniform chemical 
potentials, and so it captures qualitative, but probably not quantitative aspects of the edge-state physics (see also a 
related discussion in the end of Sec. IVip . 

When written in components, the Dirac equation corresponding to the Hamiltonian ([9]) takes the following form: 

+ + A(-) hvF{iD, + Dy)\ f^K^^ 
\ hvF {tD, - Dy) E + - A^-y Ki^K+Bj 

' E + + A(+) -hvF [iD., + Dy\ (^K.B\ 
-hvF{iD^-Dy) £: + ^(-)-A(+7 \i(jk.aJ 



0, (12) 

: 0. (13) 



Here the covariant derivative Di ~ di + {ie/hc)Ai and the shorthand notations /i^^' = jj. ± jl and A^^^ = A ± A 
were introduced (the spin index s was omitted). In each of the two sets of equations, the S-components of the wave 
function can be eliminated, 

hvF{-iD^ + Dy)ijK+A 
^^-^ = E + ^.W-Ai-) ' 

hvF {iD^ + Dy)il>K.A 
^^-^ ^ i^ + M(-)+A(+) ■ 
By taking these into account, we derive the equations for the A-components, 

{-l^Dl-PDl + l)i'K^A = 2X+^Pk^a. (16) 
{-fDl^l^Dl~l)^K^A = 2A_V/f-A, (17) 

[E + fi'-'^^Y ~ ^^^^'^^ /^O' ' = \J^(^I\^^\ is the magnetic length, and eq = ^2?iu|,|ei3|/c is the Landau 
energy scale. Note that A+ (A-) is related to the Kj^ {K-) valley. 

In the Landau gauge A = {—By, 0), the wave functions are plane waves in the x-direction. Thus, we write 

ipK+A{r,k) = -^j==u+{y,k), tpK+B{r,k) -^j==v+{y,k), (18) 
V Ztti V Ztti 

ipK^A{r, k) = —=u_{y, k), ipK^sir, k) = — ==w_(y, k). (19) 
v27r/ \2i:L 

The envelope functions u±{y, k) and v±(y, k) depend only on a single combination of the variables, ^ — y/l — kl, and 
satisfy the following equations: 

{dl~eTl + 2X±)u±{0=0, (20) 

Note that the wavevector k determines the center of the electron orbital along the y-direction, yk = kP. Then, as we 
shall see below, for a system with a ribbon geometry, e.g., Q < y < W, the condition of finite energy will be satisfied 
only for eigcnstates with wavcvectors A: in a finite range, < /c < W/l'^. This is known as the position- wavevector 
duality in the Landau gauge. (Note that the maximum value of the wavevector fc, measuring the displacement from 
either or A'_ point, is limited by the boundaries of the first Brillouin zone. However, this fact is not explicit in 
the low-energy theory.) 

The general solution to Eqs. (PO)) and (PTjl is given in terms of the parabolic cylinder (Weber) functions U (a, z) and 
V{a,z)^ 

= C, ^^"'l-^'-' u (1^, V2() + ftV (1^, V2() . (22) 
^ -CW {-i±l^. V2«) - c / + '-';'/^''V (-1±|^, V2^) . (23) 
„_,a = C,U (-i^. V2?) + C.£±i^^^F (-l±|i., V2f) , (24) 
„_,0 . c^^lilil^t, {^.V2() + C,V (i^, V-^) ^ (25) 



FIG. 1: (Color online) The lattice structure of a finite width graphene ribbon with zigzag (left panel) and armchair (right 
panel) edges. 



Note the following relations with the parabolic cylinder functions Du{z): 

U{a,z) = D_a-i/2{z), (26) 
V{a,z) = £i^±lM [sin(7ra)i?_,_i/2(z) + i?_„_i/2(-z)] . (27) 

In an infinite system, the normalizable wave functions contain only the parabolic cylinder [/-functions which are 
bound at z ^ ±oo provided a = —n — 1/2 and n is a non-negative integer. In fact, the following relation is valid: 
U{—n — 1/2, z) = 2^"/^e^^ /^i7„(z/A/2), where Hn{z) are the Hermite polynomials. In this case, the spectrum is 

determined by Aj?""^' — n with n = 0,l,2,.... 

In the case of a graphene ribbon of a finite width in the y-direction, < y < W , and with two zigzag edges parallel 
to the x-direction, the A- and i?-components of wave functions should vanish on the opposite edge o^'^i'^^ , i.e., 

^==0: u+(-fcO = w_(-fcO = 0, (28) 
y = W: v+ {W/l ~ kl) = v. [W/l - kl) = 0, (29) 

see Fig. [T] In principle, by satisfying these equations and using the wave function normalization conditions, we can 
determine all four integration constants in Eqs. ([22|l -(|25 p . For our purposes here, however, it suffices to determine 
the conditions when non-trivial, normalizable solutions exist. These will provide the dispersion spectra of all modes 
in a ribbon of graphene. The corresponding numerical analysis is presented in Sec. IIVI 

In the case of a graphene ribbon with armchair edges parallel to the y-direction, it is convenient to choose a different 
Landau gauge with {Ax,Ay) = (0,Bx). Accordingly, the solutions arc translation invariant along the j/-direction, 

i^AKir,k) - -^e^''yu+{x,k), i:BK{r,k)^^=e'^yv+{x,k), (30) 

^pAK.{r,k) = ^=e'''yu_{x,k), ^JsK-ir, k) ^ ^= e'''yv_{x, k). (31) 
\/2ttI \/2ttI 

Then, the corresponding differential equations for functions u±{x, fc), which arc rewritten in terms of the dimcnsionlcss 
variable ^ ~ x/l + fcZ, coincide with Eq. (|20|) . The expressions for the eliminated components i'±(C)i however, slightly 
differ from Eq. ()21|) . and are given by 



(32) 
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The general solutions for the w±(^) functions have the same form as in Eqs. ([22]) and ([2 

= C ^^^'l-^'-' u (1^, Voi) ,Cw{^.V-2(). (33, 
..-(0 ^ ftt/ (-i±|^, V2^) + C, ^^"'-;;^''' v V2f) , (34, 

but with ^ = x/l + kl. By using the relations in Eq. (|32p . we also obtain the explicit expression for w±(^) functions, 

. .c.a ^/2^) + .a ^ + "';'/^''' v- ^/2f) . (35, 



= ic^ ^ [/ — ^ — , V2cj + *C4V^ — ■ ^^^^ 

Since the armchair edges have lattice sites of both A and B types, the wave function should vanish both at the a; = 
and x^W lines ,^^'^'2^ 

a; = 0: u+{kl) + u^{kl) ^ v+{kl) + v^{kl) ^ 0, (37) 

x^W: u+(W/l + kl) + u^(W/l + kl) = 0, v+{W/l + kl) + v-{W/l + kl) = 0. (38) 

It is important to notice that the armchair boundary conditions mix the chiralities associated with the and if _ 
valleys, which makes the analysis more involved than in the case of zigzag edges. The details are given in Sec. |Vl 

IV. NUMERICAL RESULTS IN THE CASE OF ZIGZAG EDGES. 

Let us start from the boundary conditions at A'+ valley, see pS]) and ([29]) . They take the following explicit form: 

C ^^^'l-^'-' u (i^, -r2kl) ^ -C,V (i^, , (39, 

(-i±|^.^/2(*„ - k)l) = -c / + ''';-+'^''V (-l±|^.s/2(fo - «i) , (40, 

where fco = W/l'^ is determined by the width of the ribbon. A nontrivial solution to this set of equations exists when 



A+[/ ( i^, - V2fcA V f-i^l^, V2(fco - fc)A - U f-iil^, V2(fco - fc)A V ( -V2kl ] = 0. 



(41) 



A similar condition is derived at valley, 

>^-U i V2(fco - k)l] V f-i^, -^/2fcA - U f-i^, -^/2fcA V ( ^/2(fco - k)l] ^ 0. 



(42) 

By solving Eqs. (j4ip and (j42p numerically, we determine the dependence of dimensionless energy parameters A+ and 
A- on the wavevector k. The results for two different widths of graphene ribbons, = 5/ and W ~ 10/, are shown 
in Fig. H 

When the width of the ribbon is less than about 3 or 4 times the magnetic length I ~ 2h7K/ B\T] , we find 

that the spectra have little overlap with the usual bulk spectra, i.e., — n where n is a non- negative integer. 

Additionally, the separation between the nearest levels quickly increases with decreasing W . 

For the case W = bl, shown in the left panel of Fig. [21 only the lowest level may have a hint at the middle plateau 
developing. However, when the ribbon's width is larger than about 6 or 7 times the magnetic length, nearly flat 
plateaus are already distinguishable in the lowest levels around the central wavevector i/cp. We also find that the 
lower the level, the wider the plateau formed. 

Let us also emphasize the following special feature of the spectrum in a graphene ribbon with zigzag boundaries. As 
we see from Fig. [21 for A+ ~ (actually, E ~ — /i'+^+A^"^) and for A_ ~ (actually, E ~ — /i(~^+A(+'), dispersionless 
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surface solution o^^'^^ exist at both valleys (note that such solutions exist also in the case with no magnetic fieldf^i^). 
These solutions are bound to the /c ~ and fc ~ /cq edges for the and K- valleys, respectively. It is noticeable 
that unlike the case of a half-plane;^ they cease to be dispersionless at the opposite edges, i.e., at fc ~ fco (fc ~ 0) for 
the Kj^ {K_) valley, respectively. 

Now, by restoring the spin index, we assemble the complete spectrum of a graphene ribbon described by Hamiltonian 
(PI with dynamical order parameters proposed in Refs. I3lll32l . 




Notice that there exist eight sublevels that correspond to the lowest Landau level. Only half of these correspond 
to the hulk states, i.e., those which remain normalizable on an infinite graphene plane. The other half have wave 
functions localized only at the edges. In fact, using the properties of the parabolic cylinder functions, one can show 
that in the coordinate space, the wave functions of the additional branches of solutions (i.e., those which disappear on 
an infinite plane) are localized near either y = Q ov y = W edges of the ribbon for all values of the quantum number 
k (the non-hulk states). In other words, the earlier mentioned position- wavevector duality can be used only for the 
description of the branches of the hulk states. 

By making use of the numerical results for A-|- and A_, we can plot the actual energy spectra in the system. For the 
ribbons of widths W = 'bl and W ~ lOZ, these are presented in Figs. [3] and U) Here the choice of the order parameters 
resembles those of the solution in Eq. pTjl corresponding to the Dirac neutral point, i.e., the v ~ Q platcauj ^^i'^^ 
However, in order to lift the degeneracy of all sublevels, we also added small non-zero values for the triplet chemical 
potentials /i±. 

By considering different relative strengths of ferromagnetic and mass gaps in Figs. [3] and [H we see that there exist 
gapless edge states (whose energy vanishes at certain values of k) only when the ferromagnetic gap dominates over 

the mass gap, i.e., > |Ai^^|. (The values of wavevectors that give gapless modes are marked by the dots in 

the spectra in Fig. [31) From Fig. [21 we can see that A±(0, k) is nonnegative and approaches zero at certain values 
of the wavevector. This feature together with dispersion relations and ([33]) makes it clear that the necessary 

and sufficient condition for the existence of gapless modes is that at least one of the inequalities ■'| > |Ai^''|, 
> IaI ■'I is satisfied for at least one spin choice. 
So far, we have considered only the case with nonzero singlet Dirac masses. However, as is clear from Eqs. (|43p and 
([33]), the results for nonzero triplet Dirac masses will look exactly the same as in Figs. [31 and [31 At first sight, this might 
appear to be surprising because the symmetry properties of the two types of the masses are different: while the triplet 
masses break the SU{2)s valley symmetries, the singlet masses do not. Let us turn to the discussion of this point. 
First of all, the corresponding symmetry is exact only on an infinite plane. As to a finite width ribbon, it is explicitly 
broken by the boundary conditions there, as seen from the comparison of Eq. (|¥T|) for valley and Eq. ([1^ for A'_ 
valley. This is also obvious from the solutions for \± in Fig. [3 On the other hand, although one might expect that 
the symmetry arguments are to be approximately applicable to the states with intermediate values of the wavevector 
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FIG. 3: (Color online) Numerical results for the low-energy spectra for ribbons with zigzag edges. The ribbons' widths are 
W = 51 (left panel) and W — 101 (right panel). The ferromagnetic gaps and dynamical masses are as follows: ^± — ^O.OSeo, 
fl± — O.Oleo, A± — ±0.02eo, A± = 0. The ferromagnetic gap dominates over the mass gap, insuring the presence of gapless 
edge states (marked by dots). The electron spins of the lowest energy sublevels are marked by arrows. The spectra around the 
K+ (K-) point are shown by solid (dashed) lines. 
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FIG. 4: (Color online) Same as Fig. |3] but with the ferromagnetic gaps and dynamical masses being as follows: /i± = ^0.02eo, 
fi± = O.Oleo, A± — itO.OSeo, A± = 0. The mass gap dominates over the ferromagnetic gap, insuring the absence of gapless 
edge states. 



k in a, bulk of a ribbon, both these types of the Dirac masses still lead to the same spectra there. This puzzle is 
resolved as follows. As was already emphasized above, one should distinguish between the bulk states and non-bulk 
ones on a ribbon. The dispersion relations in the ribbon bulk for the former approximately coincide with those of the 
states on an infinite plane. In the case of the LLL, they are Esk+ = —fJ-s^^ + aI"-* and Esk^ = —fJ-s~^ + Ai'^^ 
The point is that on the ribbon, the asymmetry of this spectrum for 7^ is washed away by the inclusion of four 
additional non-bulk hranchcs, which are localized exclusively on the edges and whose energies are approximately given 
by similar expressions but with the opposite signs in front of aI^^. As for higher LLs, it is known that their energies 
do not depend on the type of Dirac masses even on an infinite plane, provided only a mass of one type is presenti^i^^ 
Thus, on a ribbon, the sublevel structure of the LLL becomes similar to that of higher LLs. 

It is instructive to compare the properties of the solutions on a finite width ribbon with those on a half-plane.— 
First of all, instead of the eight sublevels of the LLL on a ribbon, there are only six ones on a half-plane. The reason of 
that is of course connected with the fact that each edge adds two additional sublevels connected with two orientations 
of spin. Secondly, unlike the case of a ribbon, there are LLL sublevels which are dispersionless for all values of fc on a 
half-plane. This feature can be also easily understood. As the width of a ribbon goes to infinity, the edge fc = /cq 
disappears and the solutions which are bound to the edge fc = become dispersionless for all values fc > 0. 

In relation to the dispersionless modes, we would like to add the following comment. Because of level crossing, 
there is an ambiguity in the definition of a single branch of solutions in the LLL. By mixing together the two branches 
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associated with different valleys, one can construct a new branch that is approximately dispersionless for all values of 
fc, just as seen in numerical calculations!^ As should be clear from our analytical analysis, however, such a branch is 
only approximately dispersionless for intermediate values of the wavevector k. 



V. NUMERICAL RESULTS IN THE CASE OF ARMCHAIR EDGES. 

Let us now consider the case of the armchair edges. Since the armchair boundary conditions ()37|1 and (|38p mix the 
chiralities associated with the A'+ and valleys, this case is essentially more complicated for both the analytic and 
numerical analyses on a ribbon. 

The explicit form of the boundary conditions (j37p and ([38|l reads: 



^C:,U V2») + a ^ + ""l^^'*' V ^/2ti ) = 0, (45) 



C\U 



+C3 ^ + ^-^'j^- ^^"^ U [^,V2kl) + C,V (1^, V2kl) ^ 0, (46) 

+C3t/ V2ik - .0)^) + C.^±i^^l±^V^ (-i±|^, ^/2(fc - .0)^) = 0, (47) 

{-'-^,V2ik - .0)/) + C.^±J^^l-t^V (-i±|^, V2(. - k.)l) + 
^^^^B + M^^-A(+) ^ (^i^l^^ V2(fc - ko)lj + C,V {^—^. V2{k ko)l^ = 0, (48) 

where fco = —W/P. It is a homogeneous system of four linear equations for four unknown constants Ci, i = 1,2, 3, 4. 
The nontrivial solution exists when the determinant of a matrix made of the cocfhcicnts is zero. We numerically solve 
the corresponding equation to determine the spectrum of single-particle states in the graphene ribbon. 

The cases with nonzero singlet (As) and triplet (A^) Dirac masses are studied separately. Numerical results for 
singlet masses are shown in Fig. [5l We find that gapless edge states always appear when there are only Zeeman- 
like (QHF) gap and singlet Dirac masses A^, irrespective of the actual relation between their values. We conclude, 
therefore, that for the case of a singlet Dirac mass, the condition for the existence of gapless modes is less constrained 
at the armchair edges than at the zigzag ones. In fact, this result should have been expected after recalling that the 
singlet Dirac masses do not break the global valley symmetry groups SU{2)s C U{2)s. These symmetries (one for 
spin up and the other for spin down) protect the double degeneracy of the Landau levels in the bulk. The property 
of such sublevels in the LLL is that they repel in opposite directions near the edges both in the absence of singlet 
Dirac masses and when they are present. The latter leads to gapless states. (Note that in this argument we implicitly 
use the fact that the energy separation of the sublevels may become arbitrarily large. This is not true on a lattice. 
However, even on the lattice the energy separation may become much larger than the dynamical scale of the mass, 
suggesting that the conclusion is still valid in that case.) In fact, the absence of gapless modes on a ribbon with 
zigzag edges and jA^j > |/is| is a special property, which is related to the dispersionless nature of the LLL modes at 
the edges with A: ~ or fc ~ /cq- 

Let us now turn to a triplet Dirac mass. We found that the analysis of the energy spectra can be considerably 
simplified in the case of nonzero triplet Dirac masses but vanishing fls and Ag. The central observation is that in 
this case the determinant of the matrix corresponding to Eqs. (|l5l) " (P5)) can be reduced to factorized form, and the 
spectral equation becomes 



/+(fc,A)/_(fc,A)-0, 



(49) 
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FIG. 5: (Color online) Numerical results for the low-energy spectra for a ribbon with armchair edges of width W — 101 in 
the case of nonzero spin splitting and nonzero singlet masses. Nonzero dynamical parameters are as follows: ^± — ^0.02eo, 
A± — ±0.08eo (left panel) and ^± = ^O.OSeo, A± = ±0.02eo (right panel). Gapless edge states (marked by dots) are present 
in both cases. The electron spins of the lowest energy sublevels are marked by arrows. 



where, by definition, A = [{E + ^)^ — A^]/eQ and 
1 + 2A 



/±(fc,A) = U 
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± Va 
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2 ' 
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U 
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1 - 2A 



V2ik - ko)l^ V 
,V2{k - ko)l] V 



^/2fc^ j - U (-^^, V2klj V [^,V2{k - ko)l 
^ + ^^,V2kl] U ( V2kl] V f-i^, ^/2(fc - ko)l 



U - 



2 

1 + 2A 



V2{k - ko)l^ V 



2 

1 - 2A 
2 

1 + 2A 



^kl]-Ui V2kl] V ( V2{k - ko)l 



V2kl]-Ui _i±^, V2fcA V f-ii^, V2(fc - ko)l 



Restoring the spin index again, the complete energy spectrum takes the form: 



Esn{k) = -fis ± ^Jx{kl,n)€l + Al 



(50) 



(51) 



where n = 0, 1, 2, . . .. Our numerical results for several lowest branches of A versus kl are presented in Fig. [6l Note 
that X{kl, n) ~ n at intermediate values of kl. This means that the low-energy spectrum of a finite width ribbon far 
from the edges should be approximately the same as the spectrum of graphenc on an infinite plane. 

The energy spectrum in the case of triplet Dirac masses are shown in Fig. [71 As seen from the figure, the presence 
of gapless modes in the spectrum sensitively depends on the relation between the values of the Zceman-like (QHF) 
gaps and the triplet Dirac masses A^. Such modes exist only when the magnitude of the Dirac mass is less than the 
QHF gap. The physics underlying this result is clear. When |As| > \fj.s\, the valley splitting is large and there are 
two LLL branches of states with opposite signs of their energies for each direction of the spin [see Eq. (|5T|) and the 
left panel in Fig. [7]. As a result, repelling of these branches at the edges does not lead to the creation of gapless 
modes. On the other hand, in the case of |As| < \^s\, shown in the right panel of Fig. [7| [see also Eq. ([5T|) ]. the valley 
splitting is small, the energies of the two branches have the same sign, and the process of their repelling at the edges 
inevitably creates gapless modes. In essence, this is the same condition as for zigzag edges in a system with triplet 
Dirac masses studied in the previous section. 



|A^| < In, 



(52) 



VI. DISCUSSION 



The main result in this paper is establishing the criterion for the existence of gapless modes among the edge states 
in a graphenc ribbon with zigzag and armchair type edges. The method used in this paper combines analytic and 
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kl 

FIG. 6: (Color online) Numerical solutions of Eq. (|49|) for the dimensionless parameter A = [{E + — A^J/eq in the case 
of armchair edges when /i = A = 0. In the higher Landau levels, n > 1, the solutions to the equations /-(fc, A) = and 
f+{k, A) = are shown by solid and dashed lines, respectively. 
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FIG. 7: (Color online) Same as Fig. (5] but for the case of nonzero triplet masses. Nonzero dynamical parameters are as follows: 
fj,± — =p0.02eo, A± — O.OSeo (left panel) and fj,± — =p0.08eo, A± — 0.02eo (right panel). The existence of gapless modes (marked 
by dots) depends on the relative magnitude of |/i±| and |A±|. The electron spins of the lowest energy sublevels are marked by 
arrows. 



numerical approaches that allows to get a deeper insight into the nature of edge states. 

In the case of zigzag edges, gapless modes exist when the ferromagnetic (Zceman-like) gap dominates over the mass 
gap of any type, or, more formally, when any of the conditions > I Aj^'^ | are satisfied for at least one spin choice 

s = ±. This is consistent with the two limiting cases analyzed in Ref. llSl. 

For a ribbon with armchair edges, the condition for the existence of gapless modes is more involved however. In 
this case, it depends on the actual type (singlet or triplet) of the dynamical Dirac mass induced. For singlet Dirac 
masses, there are always exist gapless modes. On the other hand, for triplet masses, gapless modes exist only when 
there is a sufficiently large ferromagnetic (Zeeman-like) gap that dominates over the masses. 

One of the most interesting consequences of our finding here is a possibility of resolving the seemingly contradicting 
interpretations of the 1/ = plateau in terms of either quantum Hall metal or insulator regimesJ^ii^ As follows 
from Eq. ([TT|) , the criterion for the existence of gapless edge states at the Dirac neutral point takes the simple form 
Z + A > M (note that since a boundary of real graphene samples consists of both zigzag and armchair edges,— it 
is appropriate to use the more constrained condition for the existence of gapless modes corresponding to the zigzag 
edges). As was pointed in Refs. [siUs^ . this condition implies the existence of a critical value of the transverse 
magnetic field B^'^^\ where the insulator regime switches to the metallic one. A crude estimate for this value yields 
0.01 T < < 200 T for an effective couphng constant A taken in the interval 0.02 < A < 0.2^ As one can see, 

j^icr) ygj.y sensitive to the choice of A. In order to fix the values of A and B'-"'^ more accurately, one should utilize 
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more realistic models of graphene that incorporate consistently disorder among other thingSi ' This is a topic for 
future studies however. 

There is another possibility to explain the experimental results in Ref. [l6l As was shown in Ref. [3^ . besides the 
SI solution (jlip . there is another, triplet (T), solution around the Dirac neutral point corresponding to the 1^ = 
plateau, in the model of Refs. I3lll32l . In the T solution, while both spin up and spin down quasiparticle states have 
a triplet Dirac mass A± = M, the chemical potentials fi± are small, /i± = =pZ. Therefore, there are no gapless 
edge states for this solution and it describes the quantum Hall insulator regime. Calculating the difference of the 
free energy densities for these two solutions, it was shown^ that it is the Zeeman term which makes the 5*1 solution 
more favorable: without it, the SI and T solutions would correspond to two degenerate ground states. On the other 
hand, as was pointed out in Ref. [13 (see also Refs. [25ll47ll48l ) . there are small on-site repulsion interaction terms on 
the graphene lattice which favor the triplet solution (such terms were ignored in the model in Refs. 31,32). It would 
be interesting to figure out the role of these terms in choosing the genuine ground state in the present dynamics at 
different values of a magnetic field. 

In the future, it would be interesting to extend the present analysis by considering inhomogeneous QHF and 
MC order parameters, which should be consistently determined from the gap equation. Such inhomogeneous order 
parameters on a ribbon with zigzag edges in a magnetic field can be expected because they exist on ribbons with the 
zigzag edges in the absence of a magnetic field.— i^i^i^ The point is that as a consequence of the nearly dispersionless 
(flat) subbands, a peak in the density of states occurs near the zigzag edge, resulting in an ordered magnetic phase 
even at zero magnetic fieldi^ In the presence of a strong magnetic field, however, the MC and QHF order parameters 
alone remove the degeneracy of the dispersionless states with opposite spins (or pseudospins, related to the SU{2)s 
symmetries discussed in Sec. |II]). This means that the edge magnetism (or pseudomagnetism) may at least partially 
be captured by the homogeneous order parameters already included in this paper. Of course, it would be important 
to reexamine these arguments in more detail in future studies. 
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